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Abstract 

For m number of bosons, carrying spin (5'=1) degree of freedom, in Q number of single par- 
ticle orbitals, each triply degenerate, we introduce and analyze embedded Gaussian orthogonal 
ensemble of random matrices generated by random two-body interactions that are spin (S) scalar 
[BEG0E(2)-S'l]. The embedding algebra is C/(3) D G D Gl SO{3) with S0{3) generating spin 
S. A method for constructing the ensembles in fixed-(m, S) space has been developed. Numerical 
calculations show that the form of the fixed-(m, S) density of states is close to Gaussian and level 
fluctuations follow GOE. Propagation formulas for the fixed-(m, S) space energy centroids and 
spectral variances are derived for a general one plus two-body Hamiltonian preserving spin. In 
addition to these, we also introduce two different pairing symmetry algebras in the space defined 
by BEGOE(2)-5'l and the structure of ground states is studied for each paring symmetry. 

PACS numbers: 05.30.Jp, 05.30.-d, 03.65.Aa, 03.75.Mn, 21.60.Fw 



Corresponding author. Phone: +91 265 2434188 EXT. 211, Fax: +91 265 2423898 
E-mail address: ndchavda-apphy@msubaroda.ac.in (N.D. Chavda) 



I. INTRODUCTION 



Embedded Gaussian orthogonal ensembles of one- plus two-body interactions for finite 
isolated interacting spin-less many boson systems [denoted by BEG0E(l-|-2)] were intro- 
duced and studied in detail in the last decade in [1-5]. Going beyond spin- less boson 
systems, very recently embedded Gaussian orthogonal ensemble of random matrices for two 
spices boson systems with F-spin degree of freedom for Hamiltonians that conserve the to- 
tal F-spin of the m-boson systems [called BEG0E(l-|-2)-F] is introduced and its spectral 
properties are analyzed in detail in [6]; the F-spin for the bosons is similar to the F-spin 
in the proton-neutron interacting boson model (pnlBM) of atomic nuclei [7]. Another in- 
teresting extension of BEGOE is to a system of bosons carrying spin one {S = 1) degree of 
freedom. With random two-body interactions preserving many boson spin S then generates 
the ensemble called hereafter BEG0E(2)-S'l. In the presence of a mean- field the corre- 
sponding ensemble is BEGOE(l-|-2)-5'l. The purpose of the present paper is to introduce 
this ensemble and report results of the first analysis, both numerical and analytical, of this 
ensemble. BEGOE(l-|-2)-5'l ensemble will be useful for spinor BEG discussed in [8, 9] and 
in the analysis of IBM-3 model of atomic nuclei (here spin S is isospin T of the bosons in 
IBM-3) [10, 11]. Moreover, there is a considerable interest in analyzing a variety of embed- 
ded ensembles as they can be used as generic models for many-body chaos [12, 13] and hence 
useful to analyze BEGOE(l-|-2)-5'l ensemble. Now we will give a preview. 

In Section II, introduced is the embedded ensemble BEG0E(l+2)-S'l [also BEG0E(2)- 
5*1] for a system of m bosons in Q number of sp orbitals that are triply degenerate with total 
S'-spin being a good symmetry. A method for the numerical construction of this ensemble 
in fixed-(m, S) space is described. In Section III, embedding algebra, U{Q) [SU{S) D 
SO{3)] for BEGOE(l-|-2)-5'l is described. Section IV contains some numerical results for 
the ensemble averaged eigenvalue density, nearest neighbor spacing distribution (NNSD), 
width of the fluctuations in energy centroids and spectral variances. In addition, propagation 
formula for fixed- (m,S') energy centroids for general one- plus two-body Hamiltonians that 
preserve S and a method to propagate the spectral variances are given. In Section V, two 
types of pairing in BEGOE(l-|-2)-5'l space are introduced and some numerical results for 
ground state structure vis-a-vis the two different pairing interactions are presented. Finally, 
Section VI gives conclusions and future outlook. 
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II. DEFINITION AND CONSTRUCTION 



Let us consider a system of m (m > 2) bosons with spin 1 (5" = 1) degree of freedom and 
occupying Q number of sp levels. For convenience, in the remaining part of this section, we 
use the notation s for the spin quantum number of a single boson, s for the spin carried by 
a two boson system and for m > 2 boson systems S for the spin. Therefore, s = 1; s = 0, 1 
and 2; S = m, m — 1, . . ., 0. Similarly, the Sz ('hat' denoting operator) eigenvalue is denoted 
by TTT-s, rus and Ms respectively. Now on, the space generated by the sp levels i — 1,2, . . ., fl 
is referred as orbital space. Then the sp states of a boson are denoted by |i; s = 1, mg) with 
i — 1, 2, . . . , Q and rus — +1, and —1. With Q number of orbital degrees of freedom and 
three spin (mg) degrees of freedom, total number of sp states is A^" = 3Q. Going further, two 
boson (normalized) states that are symmetric in the total orbital x spin space are denoted 
by K^j); s, rUs) with s = lxl = 0, 1 and 2; however, for i = j only s = 0, 2 are allowed. 

For one plus two-body Hamiltonians preserving m-particle spin S, the one-body Hamil- 
tonian h{l) is defined by the sp energies e^; i — 1,2, with average spacing A, 

h{i)^Y.^,hi (1) 

i=l 

where = Ylims '^i-ms counts number of bosons in the i-th orbit. Similarly the two-body 
Hamiltonian y (2) is defined by the two-body matrix elements V^^^.^ (2) = (^{kl)s,ms \ V{2) | (ij) 
s, nis) with the two-particle spin s taking values 0, 1 and 2. These matrix elements are 
independent of the rris quantum number. The V{2) matrix in two-particle space will be 
a direct sum three matrices generated by the three ^^(2) operators respectively. Now the 
BEGOE(l+2)-51 Hamiltonian is 

{^(1 + 2)} = h{l) + Ao {V'='{2)] + Ai [V'=\2)] + A2 {v^^='(2)} (2) 

with three parameters (Aq, Ai, A2). Now, BEGOE(2)-5'l ensemble for a given (m, S) system 
is generated by defining the three parts of V{2) in two-particle space to be independent 
GOE(l)s [i.e., matrix elements are independent Gaussian variables with zero center and 
variance unity for off-diagonal matrix elements and 2 for diagonal matrix elements] and 
then propagating each member of the |if(l-|-2)| to the m-particle space with a given spin 
S by using the geometry (direct product structure) of the m-particle space. A method for 
carrying out the propagation is discussed ahead. With h{l) given by Eq. (1), the sp levels 
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will be triply degenerate with average spacing A. Without loss of generality we put A = 1 
so that the As in Eq. (2) will be in units of A. 

For generating a many-particle basis, firstly, the sp states are arranged such that the 
first Q number of sp states have rris = 1, next Q number of sp states have rris = and the 
remaining Q sp states have rris = ~1- Now, the many-particle states for m bosons can be 
obtained by distributing mi bosons in the mg = 1 sp states, m2 bosons in the mg = sp 
states and similarly, rris bosons in the mg = — 1 sp states with m — nii + m2 + m^. Thus, 
Ms — (mi — ms). Let us denote each distribution of mi bosons in mg = 1 sp states by mi, 
m2 bosons in mg = sp states by m2 and similarly, mg for ma bosons in mg = — 1 sp states. 
Configurations defined by (mi, m2, ms) will form a basis for constructing H matrix in m- 
boson space. Action of the Hamiltonian operator defined by Eq. (2) on (mi, m2, mg) basis 
states with fixed- (m, Ms — 0) generates the ensemble in (m, Ms) space. It is important to 
note that the construction of the m-particle H matrix in fixed- (m, Ms — 0) space reduces 
to the problem of BEG0E(l-|-2) for spinless boson systems and hence Eq. (4) of [1] will 
apply. For this, we need to convert the H operator into Ms representation. Two boson 
states in Ms representation can be written as |i,ms;j, mg); m^ = mg -|- mg. The the two 



particle matrix elements are V fi f2 „ ,(2) 

i,jni ;j,mi ;fc,m|^;^,ms 



i2 



It is easy to apply angular momentum algebra and derive formulas for these in terms of 
Vfjf.i{2). The final formulas are. 



^i'l;j,0;fe,l;^,o(2) 



^(l + 5i,-)(l + M 



^,0;j,0;fc,0;£,o(2) 



-i;fe,0;^,o(2) 



q L ijkl ijkl 

All other V matrix elements follow by symmetries. The fact that the sp energies e are 



independent of mg, Eq. (3) above and Eq. (4) of [1] will allow one to construct the H- 
matrix in (mi, m2, ms) basis for a given value of m and Ms — 0. Then, S'^ operator is used 
for projecting states with good S, i.e. to covert the if-matrix into direct sum of matrices 
with block matrices for each allowed S value. Matrix elements of S*^ in s = 0, 1 and 2 
space are —4, —2 and 2 respectively. This procedure has been implemented and computer 
programmes are developed. Some numerical results obtained using these programmes will 
be discussed in Section III. Let us add that the BEGOE(l+2)-5'l ensemble is defined by 
five parameters (Jl, m, Xq, Ai, A2) with As in units of A. 

III. U{n) (g) [SU{3) D 50(3)] EMBEDDING ALGEBRA 

Embedding algebra for BEG0E(l+2)-S'l is not unique and following the earlier results 
for the IBM-3 model of atomic nuclei [10, 11], it is possible to identify two algebras. They 
are: (i) [/(30) D C/(0) ® [U{3) D ^0(3)]; (ii) C/(30) D 50(30) D 50(0) O 50(3). Here 
we will consider (i) and later in Section V we will consider briefly (ii). 

Firstly, the spectrum generating algebra U{3Q) is generated by the (30)^ number of 
operators Ug{i,j) where 

(4;s=A;s=i)^ ;fc = 0,l,2 and i,i = l,2,...,0. (4) 

Note that are given in angular momentum coupled representation with k — s x s — 
0, 1, 2. Also, are one boson creation operators, b are one boson annihilation operators and 
^i;i,ms — (—1)-'^+"*^ &i;i,-ms- The quadratic Casimir invariant of U{3Q) is 

02(C/(3n)) = J]i.^(i,j)-ii'=(j,i). (5) 
Note that -U^ ^ {-lf.J{2k + l) {T'^U'^f. In terms of the number operator n, 

^=^A-XmM,m.. (6) 

we have 

C'2([/(30)) =n(n + 30- 1) . (7) 

All m-boson states transform as the symmetric irrcp {m} w.r.t. U (3f2) algebra and 

/ - \ {"^} 

{C2{U{39))\ =m(m + 3Q-l). (8) 
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Using the results given in [14] it is easy to write the generators of the algebras U (Jl) and 
SU{?,) in C/(3Q) D U{n) (g) SU{?,). The C/(Q) generators are g{i,3) where, 

= ^ (&UA;s=i)° ;i,j = l,2,...,Q (9) 

and they are in number. Similarly, SU{3) algebra is generated by the eight operators 

K = Y1 {is=A^=i)l ;fc = i,2. (10) 

i 

It is useful to mention that {hP, hg, h'^,) generate C/(3) algebra and C/(3) D SU{3). The 
quadratic Casimir invariants of U (fl) and SU (3) algebras are, 

c^mn)) = J] ^(i,j)-^(i,0, 

3'' (11) 
C2{SU{3)) ^ -Y,k'-hK 

k=l,2 

The irreps of U (Q) can be represented by Young tableaux {/} = {/i, /2, ■ ■ ■ , /f^}, J2i ft — 
However, as we are dealing with boson systems (i.e. the only allowed U (3Q) irrep being {m}), 
the irreps of U{VL) and U{3) should be represented by the same {/}. Therefore, {/} will 
be maximum of three rows. The U {Vl) and SU (3) equivalence gives a relationship between 
their quadratic Casimir invariants, 

C2{U{Sl)) = C2(f/(3)) + (i^-3)n, 

C2{U{3)) = ^'■^' = \C2{SU{3)) + ^ . ^^^^ 

fc=0,l,2 

These relations are easy to prove using Eqs. (9)- (11). Given the C/(Q) irrep {/1/2/3}, the 
corresponding SU{3) irrep in Elliott's notation [15] is given by (A//) where X — fi — f2 and 
= /2 - /3- Thus, 

{mjuisn)^ {fif2h}u{n) [{X^i)su{z)\ ; 

h + f2 + h = m, /i>/2>/3>0, (13) 

A = /i - /2, = /2 - /a ■ 
Using Eq. (13) it is easy to write, for a given m, all the allowed SU{3) and equivalently 
U{Q) irreps. Eigenvalues of C2{SU{3)) are given by 

C'2(-SC/(3)))^^''^ = C2(A/i) = [A' + //' + A/i + 3(A + //)] . (14) 



Let us add that the SU{3) algebra also has a cubic invariant C2{SU{3)) and its matrix 
elements are [16], 

C^{SU{3))y^'^^ = C^iXfi) = ^ (A - A^)(2A + + 3)(A + 2^^ + 3) . (15) 
The S0{3) subalgebra of SU{3) generates spin S. The spin generators are 

Sl^V2hl, ^ C2{SO{3)) ^ ■ S\ (S^y ^ S{S + 1) . (16) 



Given a (A /x), the allowed 5" values follow from Elliott's rules [15] and this introduces a 'i^' 
quantum number, 

K — min{\ , //), min{\ , n) —2, . . . , or 1 , 

S = max(X , yu), max{\ , /x) — 2, . . . , or 1 for X = 0, (17) 
= K,K +l,K + 2,...,K + max{X , /x) for K ^0 . 

Eq. (17) gives d(^xn)iS), the number of times a given S appears in a (A /i) irrep. Similarly the 
number of substates that belong to a U{0,) irrep {/i /2 /a} are given by dQ{fi /2 fs) where 
[17], 

dnifi) dn{fi + l) dn{fi + 2) 
da{fif2h)= dn{f2-l) dnif2) dnif2 + l) ■ (18) 
dMs - 2) daifs - 1) da{fs) 
Here, dn{{g}) = f^^^"^) and dn({«/}) = for < 0. Note that the determinant in Eq. (18) 
involves only symmetric U{Q) irreps. Using the U{3Q) D U{Q) (g) [U{3) D SO{3)] algebra, 
m bosons states can be written as |m; {/i /2 fs} a; (A /i) K S Ms); The number of a values 
is d^{fi f2 f'i)-, K values follow from Eq. (17) and —5* < Ms < S. Note that m and (A/x) 
give a unique {/i /2 /a}. Therefore i7-matrix dimension in fixed-(m, S) space is given by 

d(m,S)^ J2 Wi/2/3)rf(AM)(^), (19) 

{/i /2 /3}em 

and they will satisfy the sum rule E5(2^ + l)d{m,S) = (^"tT"^)- ^^^o, the dimen- 
sion D{m,Ms = 0) of the //-matrix in the basis discussed earlier is D{ni,Ms = 0) = 
Xlsem dim, S). For example, for {Q = 4, m = 8), the dimensions d{m, S*) for S* = — 8 are 
714, 1260, 2100, 1855, 1841, 1144, 840, 315 and 165 respectively. Similarly, for (0 = 6, m = 
10), the dimensions for 5 = - 10 are 51309, 123585, 183771, 189630, 178290, 133497, 
94347, 51645, 27027, 9009 and 3003 respectively. Because of these very large dimensions, 
numerical analysis of BEG0E(l-|-2)-Sl ensemble is quite difficult. 



IV. RESULTS FOR SPECTRAL PROPERTIES: PROPAGATION OF ENERGY 
CENTROIDS AND SPECTRAL VARIANCES 

A. Eigenvalue density and NNSD: numerical results 

Using the method described in Section II, in some examples BEG0E(2)-iS'l ensemble has 
been constructed and analyzed are eigenvalue density and spectral fluctuations. Figure 1 
presents the results for the ensemble- averaged fixed- (m, 5*) eigenvalue density p'^'^{E) for the 
BEG0E(2)-S'l ensemble defined by h{l) = in Eq. (2). We have considered a 100-member 
BEGOE(2)-5'l ensemble with m = 8 and Q, — A. The strengths of the two-body interaction 
in the s = 0, s = 1 and s = 2 channels are chosen to be equal i.e. Aq = Ai = A2. In 
the construction of the ensemble averaged eigenvalue densities, the spectra of each member 
of the ensemble are first zero centered and scaled to unit width. The eigenvalues are then 
denoted by E. Given the fixed- (mjS') eigenvalue centroids Ec{m,S) and spectral widths 
a{m, S), E — [E — Ec{m, S)\/a{m, S). Then the histograms for the density are generated by 
combining the eigenvalues E from all the members of the ensemble. In the figure, histograms 
are constructed with a bin size equal to 0.2. Results are shown in Fig. 1 for 5" = 0, 4 and 
8 values. It is clearly seen that the eigenvalue densities are close to Gaussian also the 
agreements with Edgeworth (ED) corrected Gaussians are excellent. 

The nearest neighbor spacing distribution (NNSD), which gives information about level 
repulsion, is of GOE type for spin-less BEG0E(2) [1] and BEG0E(2)-F [6]. In Fig. 2 
NNSD results are shown for BEG0E(2)-S'l with m = 8 and Q = 4 for selected spin values. 
The NNSDs are obtained by unfolding each spectrum in the ensemble, using the method 
described in [1] , with the smooth density as a corrected Gaussian with corrections involving 
up to 6th order moments of the density function. In the calculations, 80% of the eigenvalues 
(dropping 10% from both ends of the spectrum) from each member are employed. It is 
clearly seen from the figures that the NNSDs are close to the GOE (Wigner) form. 

Previously it was shown that BEG0E(l-|-2) for spinless boson systems [1, 3] and 
BEG0E(l-|-2)-F for two species boson systems [6] generate Gaussian eigenvalue densi- 
ties in the dense limit and fluctuations follow GOE in absence of the mean-field. Therefore, 
combining these with the results in Figs. 1,2, we can conclude that for finite isolated in- 
teracting boson systems the eigenvalue density will be generically of Gaussian form and 
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FIG. 1. Ensemble averaged eigenvalue density p^'^(E) vs normalized energy, 
a 100 member BEGOE(2)-5'l ensemble with 17 = 4, m = 8 and spin S=0, 4 and 8. The red cm'ves 
give Gaussian representation while the green curves are Edgeworth corrected Gaussians (ED). The 
ensemble averaged values of excess (72) parameters are as shown in figure. Note that Skewness 
71 ~ in all cases. In the plots, the state densities, for a given spin S, are normalized to dimension 
d{m,S). Note that the total dimensionality of if-matrix here is ^_5(i(m, S*) = 10234. 



fluctuations, in absence of the mean-fleld, follow GOE. As discussed in [3, 6], with mean- 
field, the interaction strength has to be larger than a critical value for the fluctuations to 
change from Poisson like to GOE. 



B. Propagation of energy centroids and spectral variances 

As the eigenvalue density is close to Gaussian, it is useful to derive formulas for energy 
centroids and spectral variances in terms of sp energies and the two-particle V{2) matrix 
elements Vf^f^i. They will also allow us to study, numerically, fluctuations in energy centroids 
and spectral variances. Simple propagation equation for the flxed-(m, S") energy centroids 
(H)"^'^ in terms of the scalars n and 5*^ operators [their eigenvalues are m and S{S + 1)] 
is not possible. This is easily seen from the fact that upto 2 bosons, we have 5 states 
{m = 0, S = 0; m = 1, S = 1; m = 2, S = 0,1,2) but only 4 scalar operators (1, n, i)?, S"^). 
For the missing operator we can use C2{SU{?>) but then only flxed-(m, (A ij)S) averages will 
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FIG. 2. Ensemble averaged Nearest Neighbor Spacing Distribution(NNSD) histogram for a 100 
member BEGOE(2)-51 with m = 8 and i7 = 4. Results are shown for the spin values S=0, 4 
and 8. Here x is in the units of local mean spacing. Results are compared with Poisson and GOE 
(Wigner) forms. 



propagate [18]. The propagation equation is 

m,(A fJ,),S 



H(l + 2] 



(m + vi2) 



m,{\,fi),S 1^ \ 1,(10),1 

= m I h{l) 



m m2 C2{\li) _S{S + l) 
" 6" ^ TS ^ 9 6 



5m 5m2 ^^(A/i) + 1) 

~6~ ^ Ts" ^ 18 ^ 6 



V{2) 



\ 2,(20),0 



2,(20),2 



(20) 



+ 



m rn? C2(A/i) 
2" ^ ~6 6 



\ 2,(01),1 



Now summing over all (A //) irreps that contain a given 5* will give ^-f^(l + 2)y . This 
is used to verify the codes we have developed for constructing BEG0E(l+2)-S'l members. 

/ ^ \ "1,5 

Propagation equation for spectral variances nif(l + 2)]^\ is more complicated. Just 

m,,(A ^),S 



as with energy centroids, it is possible to propagate the variances \ \H{1 + 2)^ 
Towards this, first it should be noted that upto m = 4, there are 19 states as shown in 
Table 1. Therefore, for propagation we need 19 5*0(3) scalars that are of maximum body 
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TABLE I. {/}, (A/ii) and S labels for m < 4 bosons and the averages of X3 and X4 operators 



m 


{/} 




(^4) 





{0} 


(00) 











1 


{1} 


(10) 


1 


5 


-25 


2 


{2} 


(20) 

















2 


21 


-147 




{11} 


(01) 


1 


-5 


-25 


3 


{3} 


(30) 


1 


9 


-81 








3 


54 


-486 




{21} 


(11) 


1 





-135 








2 





-81 




{111} 


(00) 











4 


{4} 


(40) 

















2 


33 


-363 








4 


110 


-1210 




{31} 


(21) 


1 


-7 


-121 








2 


21 


-459 








3 


18 


-246 




{22} 


(02) 

















2 


-21 


-147 




{211} 


(10) 


1 


5 


-25 



rank 4. For this the invariants n, S'^, C2{SU{3)) and C3{SU{3)) will not suffice as they will 
give only 15 scalar operators. The missing three operators can be constructed using the 
SU{3) D SO (3) integrity basis operators ^"3 and X4 that are 3— and 4— body in nature; see 
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[16, 18]. Definition of these operators are given in [16] and we call the operators given in this 
paper as X^^ and X^^{k). In the present work we have employed the following definitions, 

X,^-^Xi\ X4 = 5Xni)- (21) 
V 10 



Formulas for the averages Xj((A/i),5') = given by Eqs. (8)-(10) of [16] and 

they involve SU{3) D SO{3) reduced Wigner coefficients. Using the programmes for these, 
given in [19], averages for ^"3 and X4 in the 19 states with m < 4 are calculated and the 
results are given in Table 1. Eqs. (14) and (15) respectively will give C2(A;u) and C3(A/x). 
Propagation equation for spectral variances over fixed-(A/i.), 5* space can be written as, 

n,(A/u),5 , 

1 

Ci = 1, C2 = m, C3 = m^, C4 = m^, C5 = m^ Ce = C2(A//) , 

C7 = m C2(A A^) , Cg = C2(A //) , Cg = S{S + 1) , Cio = m S{S + 1) , (22) 

Cn = m2 S{S + 1) , C12 = S[S + 1) C2(A /x), C13 = [5(5 + 1)]^ , 

Ci4 = [C2(A/i)]2 ,Ci5 = C3(A^), C16 = mC3(A/i) , 

Ci7 = X3[(AAi),5], C18 = mX3[(A/x),5], Ci9 = X4[(AAi),5] . 

As we know {Ci) for m < 4, we can use Using ( ff^ ) for m < 4 as inputs (they 



1=1 



/ 

can be calculated by explicit construction of the Hamiltonian matrices using the method 
discussed in Section II) one can solve Eq. (22) to obtain the a^'s. Then, Eq. (22) can be 

^^m,(A^),S 

numerically X3[(A //), S\ and X4[{\ ij),S\. Their values are shown for m = 6 and 8 examples 
in Table 2. Spectral variances \H^) over fixed-S' space can be obtained easily using 



used to calculate \H'^] for any m, (A jj) and S. However we still need to evaluate 



^ \m,{\ti),S 

H ) . Let us add that there are methods [20], though much more cumbersome, 

that will give directly \H'^) ■ One such method is to use (mi, 1712, m^) configurations 
introduced in Section II and evaluate traces over these spaces. Here, trace propagation is 
simple for both H and H'^ averages and also (1711,7712,7713) configuration have a definite Ms 
value. Now, a subtraction procedure using (f/'P^("*i'"^2,TO3)^ p — 1,2 will give fixed (712, S) 
energy centroids and spectral variances. This procedure is being implemented and results 
of this will be reported elsewhere. 

Calculation of energy centroids and spectral variances for each member of the ensem- 
ble will allow us to examine the covariances in these quantities. For example, normalized 
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covariances in energy centroids is defined by 



Sn(m, S : m', S') = 1 )1 )_ ^^3) 

For (m, 5") = (m', 5") they will give information about fluctuations and in particular about 
level motion in the ensemble [1]. For {m,S) 7^ {m',S'), the covariances (cross correlations) 
are non-zero for BEGOE while they will be zero for independent GOE representation for 
the m boson Hamiltonian matrices with different m or S (with fixed Q). We have computed 
self-correlations [Eii(m, S : m, S)]^^"^ as a function of spin S for 100 member BEG0E(2)-S'l 
with {m — S,Q — A) and the results are shown in Fig. 3. It is seen that the centroid 
fluctuations are large as [En]-*^/^ ~ 28%. However, the variation with spin S is weak. We 
have also calculated in some examples the variation of the average of spectral variances with 
S and the width of the fiuctuations of the spectral widths over the ensemble. Results are 
shown in Fig. 4 for a (m = 8, = 4) system. It is clearly seen from the figure that the 
variances are almost constant for lower spins and increases for S close to the maximum 
value of S; a similar result is known for fermion systems [21]. Also, as seen from the figure, 
the width of the fluctuations in spectral widths is much smaller unlike the width of the 
fluctuations in energy centroids. Let us add that near constancy of widths is a feature of 
many-body chaos [22, 23]. 



V. PAIRING ALGEBRAS AND GROUND STATE STRUCTURE 

In the BEGOE(l-|-2)-5'l space, it is possible to identify two different pairing algebras 
(each deflning a particular type of pairing) and they follow from the results in [11, 14, 24]. 
One of them corresponds to the SO{Q) algebra in U{3n) D [U{Q) D SO{n)] ® [U{3) D 
5*0(3)] and we refer to this as SO{Q) — SU{3) pairing. The other corresponds to the 
50(30) in C/(30) D 50(30) D 50(0) 50(3). Note hat both the algebras have 50(3) 
subalgebra that generates the spin 5. Here below we will give some details of these pairing 
algebras. Inclusion of pairing Hamiltonians in BEGOE(l-|-2)-51 H will alter the structure 
of ground states and this will be discussed in Section VC. 
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BEGOE(2)-Sl: (m,Q) = (8,4), 100 members 



0.29 



0.28 



0.27 



0.26 



0.25 




FIG. 3. [Sii(m,S' : m,SY/'^ giving width of the fluctuations in energy centroids scaled to the 
spectrum width, as a function of spin S for BEG0E(2)-S'l with (m = 8, = 4). 

A. SO{a) - SU{3) pairing 

Following the results given in [11, 14, 24] it is easy to identify the Q{Q — l)/2 number of 
generators U{i,j), i < j of SO{n) in f/(3fi) D [U{n) D SO{n)] ® [U{3) D S0{3)], 



U{i,j) = \/a(JJ) [g{i,j) + a{i,j)g{j,i)] , i < j ; 



1, (^{hj)=(^U,i), (y{i,j)a{j,k) = -a{i,k). 



(24) 



Note that g{i,j) are defined in Eq. (9). The quadratic Casimir invariant of SO{Q) is. 



(25) 



Applying Eq. (24) now gives. 



i<j 

= C2{Um-Y,P,P,git,j)-git,j); 
(3i(3j = -a{i,j), for i ^ j, \l3if = 1 . 



(26) 
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BEGOE(2)-Sl: (»i,Q) = (8,4), 100 members 




0.0 0.2 



FIG. 4. Ensemble averaged variances scaled by that of maximum spin are shown for BEG0E(2)- 
51 ensembles with (m = 8,17 = 4) as a function of S/Smax- The inset figure represents ensemble 
averaged widths o"(m, S) vs S for the same example. The r.m.s. deviations (over the ensemble) in 
the widths are also shown as error bars. 



Here we have introduced /3j's and the <y{i,j) are defined in Eq. (24). Now defining the 



pairing operator "P^, = 0, 2 as 



(27) 



it is easy to see that, 



fc\t 



k=0,2;q 



C2iU{Q)) - C2{S0{Q)) - n 

-C2{SU{?,)) - C2{S0m - (fi - 4)n + - 



(28) 



In the final form above we have used Eqs. (12). Thus the pairing Hamiltonian in the 
U{?,n) D [U{9.) D SO{Q)] (g) [f/(3) D 50(3)] algebra is a sum of = and 2 pairs and it 
is simply related to the SO{VL) and S'f/(3) algebras. It is possible enumerate the irreps of 
S0{VL) given a U{VL) or equivalently S'f/(3) irrep (for a given m); see [24] and references 
therein. In terms of these irrep labels and 5'f/(3) labels (A/x), eigenvalues of H-p will follow 
from Eq. (28). This and the complimentary non-compact sp(6) pairing algebra generated by 
P^, {Vq^ , h^, /ig and h will be discussed elsewhere. For a recent review on complimentary 
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TABLE II. (A ii), S, X3 and X4 values for m = 6 and 8. 



m (A^i) S 




X4 


m (A^i) S 




X4 


6 (60) 











8 (61) 


6 


375 


-7965 


(60) 


2 


45 


-675 


(61) 


7 


340 


-5116 


(60) 


4 


150 


-2250 


(42) 











(60) 


6 


315 


-4725 


(42) 


2 





-1131 


(41) 


1 


-9 


-297 


(42) 


3 





-2046 


(41) 


2 


27 


-891 


(42) 


4 


91 


-2566 


(41) 


3 


36 


-702 


(42) 


5 


117 


-3567 


(41) 


4 


132 


-2466 


(42) 


6 


105 


-1701 


(41) 


5 


117 


-1431 


(50) 


1 


13 


-169 


(30) 


1 


9 


-81 


(50) 


3 


78 


-1014 


(30) 


3 


54 


-486 


(50) 


5 


195 


-2535 


(03) 


1 


-9 


-81 


(23) 


1 


-11 


-193 


(03) 


3 


-54 


-486 


(23) 


2 


33 


-1107 


(11) 


1 





-135 


(23) 


3 


-33 


-1314 


(11) 


2 





-81 


(23) 


4 


-44 


-1954 


(00) 











(23) 


5 


-39 


-879 


(22) 











(31) 


1 


-2 


-319 


(22) 


2 





-603 


(31) 


2 


6 


-297 


(22) 


3 





-990 


(31) 


3 


63 


-1164 


(22) 


4 





-450 


(31) 


4 


55 


-640 


8 (80) 











(04) 











(80) 


2 


57 


-1083 


(04) 


2 


-33 


-363 


(80) 


4 


190 


-3610 


(04) 


4 


-110 


-1210 


(80) 


6 


399 


-7581 


(12) 


1 


7 


-121 


(80) 


8 


684 


-12996 


(12) 


2 


-21 


-459 


(61) 


1 


-11 


-553 


(12) 


3 


-18 


-246 


(61) 


2 


33 


-1467 


(20) 











(61) 


3 


54 


-1398 


(20) 


2 


21 


-147 


(61) 


4 


166 


-3994 


(01) 


1 


-5 


-25 


(61) 


5 


171 


-2919 











Algebras see [25]. Finally, the two-particle matrix elements of H-p are V^jf — 1, V^jf — 1 
and all other matrix elements are zero. 

Before going further, it is useful to mention that the Majorana operator (M) that changes 
the space labels (i, j) in a two-particle states without changing the spin labels mg related in 
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a simple manner to C2{U{3)). Denoting the spin labels by a, (5, . . ., we have 

^ - E ^iAp {bljl^) ^ = C,{U{3)) - 3n . (29) 



B. 50(3fi) pairing 

Second pairing algebra follows from the recognition that U{3Q) admits SO{3Q) subal- 
gebra and as we will see ahead, the pairing here is generated by /c = pairs b] ■ b] alone. 
Following the results in [14] the generators of S0{3fl) are easy to identify and they are, 

u^=\i,i) ; i = 1,2,...,Q , 

V,'(i,j) = [u%j) + aii,j)i-l)'u'^{j,t)] , i<j; (30) 

|«(^,i)|^ = l, a{i,j)=a{j,i), aii, j)a{j, k) = -a{i, k) . 
The operators u'^ are defined by Eq. (4). Carrying out angular momentum algebra the 
following relation between the quadratic Casimir invariants C2{SO{30,)) and C2{U{3fl)), of 
SO{Q) and C/(3Q)), can be established using Eqs. (30) and (5), 

i i<j;k /o-^\ 

= C2{U{3n)) - ^(-l)V(z,^) •«^z,z) + X^(-l)^a(^,jV(z,j) •«'=(z,j) . 

i,k i¥'j;k 

Introducing the pairing operator P+, 

^+ - E ^+(^) = ^ E ■ ^li ' p- = (^+)' (32) 

i i 

we can prove the following relationship between C2{SO{3Q)) and the pairing Hamiltonian 
Hp = 4P+P_, 

4Hp^ AP+P_ = -n + C2{U{3n)) - C2{SO{3Q)) 

= n{n + 3Q - 2) - C2{SO{3n)) ; (33) 
7i7j = -a(i,j), for i ^ j, = 1 . 

The /3 -H- a relation is needed for the correspondence between Hp and 6*2 (5*0 (31])). Im- 
portant point now being that the three operators P+, P_ and Pq = [Q + n)/2 will form 
a SU{1, 1) algebra complimentary to 5*0(30). Thus the S0{3fl) pairing is much simpler. 
With U{3Q)) irreps being {m}, the SO{3Q) irreps are labeled by the seniority quantum 
number ou where, 

a; = m,m — 2, or 1. (34) 
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and Hp eigenvalues are 

{Hp)'^''' ^^{m-uj){m + u + 3n-2) . (35) 

The two particle matrix elements of Hp are simply Vfj^jj' — 1 and all other matrix elements 
are zero. 

C. Ground state structure 

With two different pairings in the BEG0E(l+2)-iS'l space, analysis of properties of spin 
one boson systems with the following extended Hamiltonian H^xt will be interesting and 
useful, 

{He.t} = h{l) + Xo{v^='{2)] + X,{v^=\2)] + X,{v^='{2)] ^^^^ 
+XpiH-p + Xp2Hp + Xs S'^ . 
As an example, shown in Fig. 5 are the results for the probability that the ground state spin 
is 5" = Smax = m. It is seen that even with strong random interaction the probability is not 
100% and with increasing Xs, the probability rapidly comes down to zero. The situation 
here is different from the result seen for BEG0E(l+2)-F where strong enough interaction 
generates states with maximum spin with 100% probability. With S0{3fl) pairing the drop 
in probability is more rapid compared to the situation with 5'0(Q) — SU{3) pairing. In 
addition, we have also calculated the expectation values of the two pairing Hamiltonians 
and C2{SU{3)) and the results are shown in Fig. 6. We have considered BEG0E(l+2) 
Hamiltonian defined by Eq. (2) with Aq = Ai = A2 = A = 0.2, i.e. in the region of chaos 
generated by random two-body interactions in the presence of a mean-field. It is seen that 
the expectation values are largest near the ground states and then decrease as we move 
towards the center of the spectrum. Due to finiteness of the model space, the curves are 
essentially symmetric around the center. The calculated results are in good agreement with 
the prediction [6] that for boson systems (just as it was well verified for fermion systems), 
expectation values will be ratio of Gaussiaus; see Section VI and Eq. (43) in [6]. Results 
in the figure show that with repulsive pairing, ground states will be dominated by low 
seniority structure. Also, with random interactions, there is no clear distinction between 
the two different pairing structures. Variation with the A parameter for larger systems (Q 
and m large) may show the difference but these calculations are not attempted as the matrix 
dimensions will be very large. This exercise will be attempted in future. 
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FIG. 5. (a) Probability for ground state to have maximum spin as a function of A5 for a 250 member 
(m = 6,17 = 4) BEGOE(l+2)-51 system with H defined by Eq. (36). In ah the calculations 
Aq = Ai = A2 = A and the results are shown for A = 0.2 and 0.5. (a) results with SO{Q) — SU{3) 
pairing parameter Api = 0.3 and SO{30,) pairing parameter Ap2 = in Eq. (36). (b) Same as (a) 
but with SO{Q) — SU{3) pairing parameter Api = and S0{3^}) pairing parameter Ap2 = 0.3 in 
Eq. (36). 



VI. CONCLUSIONS AND FUTURE OUTLOOK 



Introduced in this paper is the embedded Gaussian orthogonal ensemble of random ma- 
trices generated by random two-body interactions in presence of a mean-field for spin one 
boson systems. Presented are some first analytical and numerical results for this ensemble. 
Due to large fixed-(m, = 0) matrix dimensions, only restricted numerical calculations 
(with dimensions less than 10000) could be carried out at present. Some results for spectral 
properties including the form of eigenvalue density close to Gaussian, NNSD following GOE 
for sufficiently strong interaction strength and also for lowest two moments of the two point 
function are presented. It is possible to deal with much larger space if we use direct con- 
struction of H matrix in a good 5* basis. This is being attempted and using this in future a 
more detailed study with much larger size examples will be reported. Preliminary aspects 
of one of the embedding algebras SU (Q) ® SU (3) and also two pairing algebras in the space 
defining BEG0E(H-2)-S'l are discussed in the paper. More detailed study of the effects of 
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Q=4, m=6 : ^=0.2 

T— I— I— I— I— I— I— I— I— I— I— i-q L I I I I I I I I I I I I 




-4 -2 2 4 -4 -2 2 4 

A 

FIG. 6. Expectation values of the two pairing Hamiltonians and C2{SU{'i)) vs E for a 250 member 
BEG0E(l+2) systems with H defined by Eq. (2) and (fi = 4, m = 6). Results are shown for spins 
5 = and 5 = 4. (a) expectation values of H-p, (b) expectation values of Hp and (c) expectation 
value of C2{SU {?>)). Ensemble averaged results are shown by histograms while continuous curves 
are ratio of Gaussians given by EGOE theory [6]. See text for further details. 

random interactions in presence of the two pairing interactions will be discussed elsewhere. 
Extension of BEGOE(2)-51 to BEGUE(2)-51 and to the more restricted BEGUE(2)-Sf/(3) 
with H preserving SU{'i) symmetry for spin one boson systems are possible; see [26] for pre- 
liminary results for BEGUE(2)-S't/(3). Finally, applications of BEG0E(l+2)-S'l ensemble 
to spin one BEG should be possible in future. 
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